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SUMMARY 

By use of the assumption that the pressure-volume relation is 
linear, a solution to the problem of designing a cascade for a 
given turning and with a prescribed velocity distribution along 
the blade in a potential flow of a compressible perfect fluid was 
obtained by a method of correspondence between potential flows of 
compressible and incompressible fluids. 

If the prescribed velocity distribution is not theoretically 
attainable, the method gives a way of modifying the distribution so 
as to obtain a physically significant blade shape. 


INTRODUCTION 

In order to control boundary- layer growth, transition, and 
separation in the design of a cascade for a given turning, it is 
advantageous to prescribe the velocity along the blade as a func- 
tion of the arc length along the blade and then to compute the 
blade shape. For the case of an incompressible fluid, several 
solutions to this problem have been obtained. (See references 1 

to 3.) 

7 « 

A similar solution for the two-dimensional potential flow of 
a compressible perfect fluid has been developed at the NACA Lewis 
laboratory. This solution is based on the assumption that the 
pressure-volume relation is given by a linear approximation to the 
isentropic curve instead of the true curve. The flow pattern of 
the compressible fluid is obtained by a transformation from a cor- 
responding flow of an incompressible fluid using the transformation 
developed by Lin (reference 4). 

The method of solution consists in using the free-stream 
velocities upstream and downstream of the cascade and the prescribed 
dimensionless velocity distribution along the blade to select a 
suitable incompressible potential flow about the unit circle and 
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then to determine the mapping function that transforms this incom- 
pressible flow into a compressible flow about a cascade of airfoils. 
The image of the unit circle under this mapping gives the cascade 
with the prescribed velocity distribution along the blade, provided 
the velocity distribution is theoretically attainable. If the veloc- 
ity distribution is unattainable, methods are given for modifying the 
distribution so that a physically significant profile is obtained. 


SYMBOLS 

The following symbols are used in this report: 


A,B,C^,C2>D 

complex constants 

a l> a 2 

location of complex sources in £ -plane 

c(t) 

function of £ defined by equation (44) 

d 

spacing of cascade 

F(£) 

complex potential function (incompressible flow) 

tit) 

regular function of £ 

g( £) 

regular function of £ 

HU) 

regular function of £ 

Be H(£) 

function of £ defined by equation (59) 

Im 

imaginary part 

K(0) 

function of 0 defined by equation (49) 

t 

K 1 

2<li 

constant equal to 

1 + f\j 1 + q^ 

K 2 

2q 2 

constant equal to 

1 + V 1 + ^2 

k 

constant defined by equation (33) 

n 

number determined by included trailing-edge angle 
of blade 
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P 

Q(s) 

q 




V 


ia l 

ia 2 


pressure 

auxiliary function of s 

magnitude of dimensionless velocity in compressible- 
flow plane (ratio of actual velocity to stagnation 
velocity of sound) 

dimensionless velocity upstream of cascade 


. dimensionless velocity downstream of cascade 


E 

Be 

r 

s 

v(e) 

z = x + iy 
a 

r 

7 

5 

£ « I + ii) 
e 

p 

T 

CP 

M> 


region in £ -plane defined by |£| > 1 
real part 

number defined by equation (23) 
arc length along blade 

velocity on unit circle (incompressible flow) 

complex variable (compressible -flow plane) 

angle of velocity in compressible flow (measured from 
positive x-axis) 

circulation (positive counterclockwise) 
ratio of specific heats 
included trailing- edge angle of blade 
complex variable (incompressible-flow plane) 
circle angle (incompressible- flow plane) 
auxiliary variable defined by equation (37) 
density 

variable of integration 
velocity potential 
stream function 
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Subscripts : 

c compressible flow 

i incompressible flow 

n leading edge 

t trailing edge 

Prime indicates a derivative. 


THEORY OF METHOD 

In reference 4, Lin has shown that if the pressure- density- 
relation is 



( 1 ) 


then the compressible potential flow about a cascade of blades can 
be obtained by transforming the incompressible flow about the uni t 
circle in the following manner. The complex potential function F(£) 
for the incompressible flow due to two complex sources at £ = a-i 
and £ = a 2 outside the unit circle. |£| = 1 is 1 


F(£) = -A log e (£-a]J + A log e (£- + B log e (£-a 2 ) + B log e D 

( 2 ) 

where A and B are complex constants with Be A >0 and 
Be A = - Be B, and D is an arbitrary complex constant. The bar 
indicates the complex conjugate. The mapping between the z -plane 
and the £ -plane defined by 


dz = g(£) a-ai)' 1 (t-ag)" 1 it - |[f’U)] 2 [gtt)]" 1 (t-e^) (J-ag) it 

(3) 

gives a compressible flow with the linear pres sure- volume relation 
past a straight cascade of identical blades in the z-plane with 
the velocity potential op c and stream function \|/ given by 

<P C + 1 % = F(0 


(4) 


28TI 
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provided that g( £) is chosen to satisfy the following require- 
ments: 

(a) The function g(£) is regular in closed region R ^ 
defined by |£| si. 


(b) The function g(£) ^ 0 in R, except possibly at 
one point on the circle where F'(£) =0. (The order of 
the zero not to exceed 1 . ) 


(c) Along the circle 


|tl - 1, <f> 

l . 


> (5) 


dz = 0. 


(d) The function g(£) satisfies the inequality 
-1 


[F'(£)].[g(£)] 


U-*l) (t-&2) 


<2 in R. 


J 


The magnitude q and the direction a of the dimensionless 
velocity at any point in the z-plane are given by 


2q 


-ia 


+V 1 + 


iTD 


( 6 ) 


In order to use this transformation in designing a blade with 
a prescribed dimensionless velocity distribution along the blade 
in a cascade, the prescribed conditions are used to select a suit- 
able incompressible flow about the unit circle and to determine the 
function g(£). . 


The prescribed conditions are the velocity distribution on the 

ia-i 

airfoil, the upstream velocity q 1 e and the downstream velocity 
q 2 e ia2 . The upstream and downstream velocities are related by the 
isentropic-flow equations with y = -1. This relation is 


q 2 2 cos 2 a 2 


1+q. 


2 2 
q 1 cos 

I 2 

l+q x 


(7) 


where the axis of the cascade has been taken along the y-axis for 
convenience and the flow is from left to right. (See fig. 1.) 
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Flow in Circle Plane 

The flow of an incompressible fluid about the unit circle is 
selected by determining the constants A, B, a^, and ag in 
the complex potential 

F(£) = A log e (t-a 1 ) + A log e + B log e (t-a 2 ) + B log e + D 

( 8 ) 

from the given conditions . The constants A and B are obtained 
from the upstream and downstream velocities and the circulation 
and then a^ and a^ are determined by the range of the potential 

on the airfoil. 


Circulation and cascade spacing. - The magnitude of the pre- 
scribed dimensionless velocity q along the airfoil is given as a 
function of the arc lengths [q = q(s)] where the total arc length 
is taken to be 2it and is measured from the trailing edge along the 
lower surface. If Q(s) is defined by 


Q(s) = - q(s) 0 < s < s n 

Q(s) = q(s) s n < s < 2it 


where s Q is the leading-edge stagnation point, then 


cp c (s) = / Q(s) ds 

do 


r 

c 



(ID 


The circulation and the spacing of the cascade are related by 



q^ sin s * n a 2 


( 12 ) 


where d is the spacing. The quantities T c , q 1 , q g , and 

ag are known so that the spacing is determined. 


1183 
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Determination of A. - The value of d from equation (12) is 
used to evaluate A and B because the spacing is also given by 
the absolute value of £ dz taken along a path around a^ or ag . 

(See fig. 2.) The axis of the cascade has been taken along the 
y-axis so that 



The second equality comes from the fact that the residues at infin- 
ity of 

siiHK)' 1 u-a 2 r l 

and 

% 

in the expression for dz in equation (3) are zero and consequently 



where c is the unit circle. But by equation 5(c) 


so that 




(14) 


(15) 


The evaluation of equation (13) in terms of the potential F(£) 
is 
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But 


F'(£) = 


B 


B 


- + r- + -f + — 

i i - ~ a 2 t . JL 


t - ^ 
a 2 


(17) 


so that equation (16) reduces to 


id = 2ni 


g(a 1 ) 1 

( a l-a 2 ) ” 4 


A (a l-°2> 
g(a x ) 


2ni 


(18) 


= 2rti 


g(a x ) j2 (a^ag) 




L (a l-*2> 4 g (*!> J 


At £ = a^, equation (6) becomes 


2 ^1 “iOi A(a 1 -a 2 ) 


1 + 




e 


g(a n ) 


(19) 


which on writing 


2q, 


Kt = 


1 = p=^ 

1 + Ajl+q^ 


reduces to 


g(a 1 ) A 


( a l-a 2 ) K 1 


( 20 ) 


gives 


Substitution of the values from equation (20) in equation (19) 


id = 2ni 


l(± e lai +S. e la l 


i K- 


( 21 ) 


Hence, the bracketed expression in equation (21) must be a real 
number and 

4A+AK] 2 _ -icu 
4K, = r ® 


( 22 ) 
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where 


2 _ 


r = 


'4+Ki 2 \ 2 /4-Kj. 2 f 

Re A ) + I —rz . — Im A ) 


, 4K 1 


4K, 


(23) 


From equation (22) 


(4-K-, ) Im A 

= - tan a]_ 

(4+K^) Re A 


or 


4+K x 


Im A = - 77 Ee A tan 


4-K-. 


°1 


Im A = 




Ee A tan 


°1 


(24) 

(25) 


which gives the relation between Ee A and Im A. 

Substitution of the value of Im A from equation (24) in 
equation (23) yields 


2 2 

/4+K-, 2 \ /4+ K x 2 \ 2 2 

r =l-iKrJ Ee Ee Atan °1 


4+K l)p A 2 

rsrv Re A sec °i 


or 


r = 


4-hKj 


Re A 


sec 


“1 


Hence 


/ 4+K l , , N 

d = 2 jt I — Re A I sec a^| 


( 26 ) 
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Substitution of the value of d from equation (26) in equa- 
tion (12) gives 


2 n 




F o 

q x sin ct^ - q 2 sin 


or 


Ee A 


4K^ r c cos OqJ 

p 

2n(4+K^ )(q^ sin - q 2 sin ctg) 


(27) 


and Re A is now determined. By use of this value of Ee A in 
equation (25), Im A is obtained. Hence A is completely given 
by equations (27) and (25). 


Determination of B. - From equation (13) 


id 


dz 






id = - 2iri 


g(a 2 ) x b 2 (a 2 -a x ) 2jti 


+ v 


( a 2“ a l) 4 g ^ a 2^ 


= - 2iti 


Be 


ia 2 2rt 


Kc 


i _ ia? 
4 BK 2 e ^ 




where 
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The bracketed expression must be real, so that 

(4-K 2 2 ) Im B 

= - tan ao 

(4+K 2 2 ) Be B 

But 

Be B = - Be A 

and equation (28) can be written 

4+K 2 2 

Im B = Be A tan cu> 


(28) 


(29) 



Be A tan 


“2 


(30) 


Consequently, B is determined by equations (29) and (30) because 
Be A is known from equation (27). 


Determination of a^ and a 2 . - After A and B are known, 
the points a^ and a 2 are to be selected to satisfy the single 

condition that the range of potential on the circle must equal the 
range of potential on the airfoil, that is. 


% < 2 *) - %( s n ) = F f i(e t +2 *)] 


F(e iSn ) 


(31) 


where 0^ and 9 n are the trailing-edge and leading- edge stagna- 
tion angles, respectively. This condition is only one equation in 
two unknowns, aj_ and a 2 ; consequently, the values of a^ and 

a 2 are not uniquely determined. By imposing an additional restric 
tion that a^ and a 2 are real and 

a i = - a 2 (32) 

unique values are obtained for a^ and a 2 in all cases. In 

particular problems, however, some other restriction may be more 
useful, such as assigning a definite value for a^ and computing 
a 2 . 
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With the restriction given in equation (32), it is possible 
to express 0^ and 0 n in terms of a^ and substitute these 

values in equation (31), but the resulting equation cannot be 
. solved explicitly for a^. 


One method for obtaining a^ is as follows: Let 


a l - 



(33) 


where k> 0. Then equation (8) becomes 

F(|;) = A log e (£-e k ) + A log e (£-e" k ) +B log e (£+e k ) + B log e (£-e" k ) . + D 


or, for points on the unit circle 


t- 


10 


e 


P(e 10 ) = cpjfe) + 1^(0) = Be A log e 


(e le -e k )(e ie -e- fc ) 

(e ie + e k )(e 16 + e- k ) 


+ 


i0 k i0 k 

i Im A log e e i0 ' e _ - k + i Im B log e + D (34) 

Hence, cp(0) may be written in the form 

cp(0) = -2 He A tanh"'*' COS , + (im A + Im B) tan ^ ft - + 

■ cosh k tanh k 

1 Q-fn A 1 COS 

(Im A - Im B) tan' 1 + 2 Be A tartT 1 - 

(L. A + a. B) tan' 1 - (1« A - Im B) tan 

(35) 


cp i (0 t ) = 0 

sin 0 < n 

sinh k 2 


where D has been chosen to make 
vention is 


and the angle con- 
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and tan"^ is taken in the same quadrant and same direction 

tanh k 

as 0. 

The velocity on the circle v(0) is 

_ /n \ 2Ee A sin 0 cosh k . (im A + Im B) tanh k sec 2 0 

“ 5 o + 5 — 5 + 

tanh 6 k + tan 6 0 


2 2 
cosh k - cos 0 


(Im A - Im B) 


slnh k cos 0 
sinh 2 k + sin 2 0 


= t „ sttt 2Be A sin 0 cosh k + 

(cosh 2k - cos 20) L 

(im A-Im B) cos 0 sinh k + (Im A + Im B) sinh k cosh kj 
Equation (36) can be further simplified by defining X as 


(36) 


tan X = 


_ (im A - Im B) sinh k 


2Be A cosh k 


‘I <X< f 


Then 

v(0) = 


_ 2 x/4Re 2 A cosh 2 k + (im A-Im B) 2 sinh 2 k 


cosh 2k - cos 20 


sin 0 cos X + 


cos 0 sin X + 


(Im A + Im B) sinh k cosh k 


/j4Re 2 A cosh 2 k + (im A - Im B)^ sinh^lcj 

sin (0 + X) + 


2 ^4Re 2 A cosh 2 k + (Im A - Im B) 2 sinh 2 k 


cosh 2k - cos 20 


(im A ± Im B) cosh k sinh k 

^4 Re 2 A cosh 2 k + (im A-Im B) 2 sinh 2_ k 


(38) 


The stagnation angles 0^. and 0 Q are therefore the roots of the 
equation 
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sin(e+A) = -(Im An-ImB) sinh k cosh k (39) 

/\fl Re^A cosh^ k + (im A - Im B)^ sinh^*"k 

The desired value of k is obtained as follows: 

(1) Assume a value of k. 

(2) Compute X by equation (37). 

(3) Obtain 0^. and 0 Q from equation (39). 

(4) Compute ^(0^+2^) - ^(Sq)* 

(5) Repeat (l) to (4) several times to obtain a plot of 
cp i (0 t +2w) - qp^(0 n ) as a function of k.. 

(6) Interpolate to obtain k such that 

<*^(0^2*) - cp i (0 n ) = <p c (2n) - <?P c (8 n ) 


With k determined the flow about the circle is known. The 
potential cp^(0) and velocity v(0) for points on the circle are 
given by equations (35) and (38), respectively. 


Function g(£) 

The function g(£ ) can be computed for points on the unit 
circle by using the prescribed velocity on the airfoil and the 
velocity on the unit circle to determine the real part of g( £) . 
The imaginary part of g( 0 can then be computed by Poisson' s 
integral. Because of the restrictions imposed by the given con- 
ditions, however, g(£) is actually obtained in a slightly dif- 
ferent maimer, as shown in the following sections. 


Airfoil with pointed trailing edge. - If an airfoil with a 
pointed trailing edge is desired, then g(£) must vanish at the 

trailing-edge stagnation point £ = e i0t . Hence, g( £) can be 
written in the form 


g(£) = 1 




(40) 
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where f ( £) is regular in the exterior of the unit circle and 



(41) 


where 6 is the included trailing-edge angle of the airfoil. 
(See reference 4.) 


Values of g(±e k ). - Because the velocities are given for the 
compressible flow upstream and downstream of the cascade, the 
value of g(£) at £ = ±e^ is determined from equation (6), 


g (e*) - 


"\ 


ia„ 




, k, - 2e k Be"'2 

g( ' e) ~ y 


(42) 


In order that g(£) have these values, f(£) is written in 
the form 


where 


■ e«) + H(t) 

■i z 


c(£) - I (i + j) log. 


if 1 - 7” V°S. 


2Ae 


ict^ + k 


r r°* 


[ K , ( X - e igt "^)0 
-2Be la S + k 


4 i9 t -k^n 


(43) 


Kg (l+e 

and H( £) is regular in the exterior of the unit circle with 


(44) 


lim i H(£) = 0 (45) 

The restriction on H(£) imposed by equation (45) is neces- 
sary so that f(£) (equation (43)) will be regular. By use of 
equation (43), g(£) is expressed as 
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g(t) = 1 - 


i0 t \ 


cd) 


(e 2k -t 2 )(t 2 -e~ 2k ) 

i 2 


H(D (46) 


and g(£) will be known when H( £) is determined. For the actual 
computation of the blade shape, only the values of g(£) on the 
unit circle are needed-. Hence, it is only necessary to compute H (0 
for points on the circle. If desired, the values of H (£) for any 
point in the exterior of the circle can be obtained from the values 
on the circle by Poisson's integral. 


Determination of Be H on the circle. - By equation (4), the 
potentials cp c (s) and cp^©) are equal at. corresponding points. 


Thus, by matching these potentials a correspondence is established 
between points along the airfoil &rc and the circle angles; that is, 
s = s(0). By use of this correspondence, the magnitude of the pre- 
scribed velocity along the airfoil is obtained as a function of the 
circle angle q = q(0). Hence, by taking absolute values of equa- 
tion (6) 


2q(0) 

i + Vi+q(e) 2 


F’(e i0 )l 

f ie k 

[e -e 

w i0 k» 
)(e +e ) 


g(e l0 )| 



(47) 


for points on the circle. Substitution of the value of g(£ ) from 
equation (46) with £ = e 1 ® and replacing F'(e i0 ) by the velocity 
v(0) (equation (38)) on the circle give 


2q(6) = 

1 + A/l+q(0) 2 ' i 

e k |v(0)| (2 cosh 2k - 2 cos 20) 2 

[2 - 2 cos (9 t -e)p el? e C(ei6) + < 2 coeh 2k - 2 cos . 2e > Ee H(eie l 


or, with the equation solved for Be H(0), 

| v(0 ) | (2 cosh 2k-2 cos 20) 


log. 


Be H(e i0 ) * 


K(0) [ 2-2 cos (0 t - 0)J 2 J 


- Be C(e i0 ) + k 


2 cosh 2k - 2 cos 20 


< 


( 48 ) 
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where 


K(e) = 


1 + 


2q(0) 

*] l+q(0) 5 


(49) 


Restrictions on Be H(e* 0 ). - Equation (45) imposes restric- 
tions on the values of Be H(e i0 ), as shown hy writing H(£) in 
the form 


h(U + J + p + • • • 

For points on the circle, equation (50) becomes 


(50) 


H(e i0 ) = Be H(e 1P ) + i Im H(e i ' 7 ) 

=s Re hQ + E (Be h, cos J0 + Im hj sin J0) + 


i0i 


.i0. 


Im hQ + y ] (Im hj cos J0 - Be hj sin J0 


(51) 


Equation (51) is a Fourier expansion and 

P2rt 

Be hn = - / Be H(e i0 ) d 0 

*Jo 


(52) 


■»2rt 


Be h i - ^ f Re H(e ie ) cos 0 d0 


(53) 


*2n 

Im hi = ~ / Be H(e 10 ) sin 0 d0 
**" 2n 


(54) 


But equation (45) requires that 

Be h Q = Im h Q = Be h^ = Im h^ = 0 


(55) 
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Consequently, 


Re H(e i0 ) 


must satisfy the equations 



H(e i0 )d0 = 0 



H(e* 0 )cos 0 d0 = 0 



H(e i0 )sin 0 d0 = 0 


(56) 

(57) 

(58) 


Adjustment of Re H(e^ 0 ) . - If the values of Re H(e^ 0 ) from 
equation (48) do not satisfy equations (56), (57), and (58), the 
values must be adjusted until the conditions are satisfied. One 
method for adjusting the function is to define Re H(e^) by 


Re H(e i0 ) = Re H(e i0 ) - | Re h Q - 2 Re h x cos 0 - 2 Im h L sin 0 

(59) 


where Re hQ, Re h^, and Im h-j_ are given by equations (52), 

(53), and (54), respectively. The modified function Re E(e* 0 ) 
will then satisfy equations (56), (57), and (58). This method of 
modifying Re H(e^®), however, changes the velocity distribution 
all along the profile and, if the correction terms in equation (59) 
are not small, these changes in the velocity may be extensive 
because 




+ 2 Re h^ cos 0 + 2 Im h^ sin 6j 


where q denotes the new velocity. In some cases, consequently, 

Re H(e^ 0 ) can best be adjusted to satisfy the requirements by 
adding to Re H(e i0 ) odd and even functions that have nonzero 
values only in small neighborhoods of the points 0=0 and 
0 * - it. The particular functions to be added to Re H(e* 0 ) and 
their range of values depends on the specific problem; no general 
method can be given for determining the functions. 
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Determination of Im H(e i0 ). - After Ee H(e i0 ) satisfying 

equations (56), (57), and (58) is obtained, the function Im H(e i0 ) 
is given by Poisson's integral (reference 5) 


Im H(e 10 ) = \- 



cot 


(¥) 


dT 


(60) 


where the constant term in Poisson's integral has been taken as 
zero so that 


2 n 

•Im H(e i0 ) d0 = 0 

as required by equation (55). Hence, H(e i0 ) is determined for 
points on the unit circle by 

•• H(e i0 ) = Ee H(e i0 ) + i Im H(e i0 ) 

Adjustment of g(e* 0 ). - By use of these values of H(e* 0 ) 
in equation (46), g(£) is determined for points on the circle, 

g(e i6 ) - [l e [c(e i0 ) + ( 2 cosh 2k - 2 cos 20) H(e 10 )] 

(61) 

Because of the adjustments in- H(e^ 0 ), g(e^ 0 ) may no longer 

satisfy condition (5d) . If g(e^ 0 ) does not satisfy the inequality 
for points on the circle, then the values of g(£) can be adjusted 
to satisfy the inequality by changing the second or higher harmonic 
terms in H( £ ) or by other methods . It should be noted that if 
the velocity along the profile is finite then g(£) satisfies the 
inequality. In fact, if the prescribed conditions are theoretically 
attainable, then no modification is necessary, not even in Ee H(e^). 



Blade Coordinates 

By use of the values of g(e i0 ) that satisfy all conditions, 
the blade coordinates are obtained from integration of equation (3), 
that is 
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g(t)(^f dt - i J - (J) ^ (|) 1 ( I + - )d J 


which on replacing F'(e*®) by v(0)e _: *-(® + 2) and writing 


g(. 18 )-= gl («) e lg2(S) 


reduces to 


2 = 


(e 


2i0 


• e 2 *)- 1 


_jief 

4g x (0) 


(e 2i0 -e 2k ) 


{® + | + S2( 0 )] 


d0 


(62) 


COMPUTATIONAL PROCEDURE 

An outline of the procedure for computing the blade shape is 
as follows: 

(1) Obtain cp c (s) and T c from equations (10) and (11), 
respectively. 

(2) Compute Ee A, Im A, and Im B by equations (27), (25), 
and (30), respectively. 

(3) Obtain k as outlined in the text. Compute q^(0) and 
v(0) by equations (35) and (38), respectively. 

(4) Plot cp c ( 8 ) and q>^(0). By comparing the abscissas for 

equal values of these potentials, obtain s as a function of 0, 
which permits writing the prescribed velocity q as a function of 

0, q « q(0). 

(5) Compute Ee H(e^) by equation (48) and determine Ee hQ, 
Ee h^, and Im h^ by equations (52), (53), and (54), respectively. 
If these values are not zero, then adjust Re H(e i ®) either by 
equation (59) or by addition of functions so that Re H(e*®) satis- 
fies equations (56), (57), and (58). 

. (6) Obtain Im H(e^) by equation (60) using the adjusted 
values of Ee H(e ie ). 
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(7) Obtain g(e*®) by equation (61). The function g(e i0 ) 
must satisfy inequality (5d) for points on the circle. If gCe 1 ®) 

does not satisfy the inequality, adjust g(e i0 ), as suggested in 
the text. 

(8) After g(e ie ) has been adjusted to satisfy all conditions, 
the blade shape is obtained by integrating equation (62). 


DISCUSSION 

The magnitude of the dimensionless velocity along the blade 
cannot be entirely prescribed arbitrarily as a function of the arc 
length, but is subject to some restrictions in addition to the con- 
ditions imposed on H(£) previously discussed. The magnitude 
must be finite everywhere along the profile and by the method given 
here the velocity can be zero in at most two places - the leading- 
edge and trailing-edge stagnation points. By a limiting process, 
however, the method can be extended to provide for additional stag- 
nation points. The zero of q at the trailing edge is of the 

s 

order — , where 8 is the included trailing-edge angle of the 

Jt 

blade. Thus, for a cusp at the tail, 8 is zero and q need not 
be zero at the trailing edge. 

Another restriction is imposed on the velocity distribution 
when the spacing of the cascade, as well as the turning, is speci- 
fied in advance because the distribution must be selected so that 
r c will satisfy equation ( 12 ) . 

If a velocity distribution is selected to satisfy these con- 
ditions but otherwise is arbitrary, the resulting profile may not 
' be a physically real blade but may result in a blade with zero or 
negative thickness in some portions of the blade. The negative 
or zero thickness is caused by specifying too low velocities along 
parts of the blade and a physically real blade can be obtained by 
increasing the prescribed velocity along the blade. 


CONCLUSION 

By use of the assumption that the pressure-volume relation 
is linear, a method has been given for computing the blade shape 
in a straight cascade of identical blades having a prescribed 
velocity distribution along the blade and given upstream and 
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downstream velocities in a potential flow of a compressible per- 
fect fluid. If the prescribed velocity is not theoretically 
realizable, the method gives a way of modifying the distribution 
so as to obtain a blade shape. Whether the resulting blade is 
practical will depend on other considerations. The applicability 
of the method is limited only by the accuracy of the linear approxi- 
mation to the pressure-volume relation. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 7, 1949. 
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Figure 2 ♦ - Paths of integration in C-plane. 
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